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We discuss critical gravitational collapse on the threshold of apparent horizon formation as a 
model both for the discussion of global aspects of critical collapse and for numerical studies in 
a compactified context. For our matter model we choose a self-gravitating massless scalar field 
in spherical symmetry, which has been studied extensively in the critical collapse literature. Our 
evolution system is based on Bondi coordinates, the mass function is used as an evolution variable 
to ensure regularity at null infinity. We compute radiation quantities like the Bondi mass and news 
function and find that they reflect the discretely self-similar (DSS) behavior. Surprisingly, the period 
of radiation at null infinity is related to the formal result for the leading quasi-normal mode of a 
black hole with rapidly decreasing mass. Furthermore, our investigations shed some light on global 
versus local issues in critical collapse, and the validity and usefulness of the concept of null infinity 
when predicting detector signals. 

PACS numbers: 04.25.Dm, 04.20.Ha, 04.20.Dw, 04.30.-w 

I. INTRODUCTION 

In this paper we present a numerical study of scalar field critical collapse on the threshold of singularity formation. 
We consider Einstein's equations with a minimally coupled massless scalar field (j): 

Gab = 87rr,b = Va(/'Vfc(/)-igafcV,(/)V^(/., (1.1) 

Og^ - 0, (1.2) 

We restrict ourselves to the case of spherical symmetry, which has been studied extensively with numerical and 
analytical techniques since the discovery of critical collapse by Choptuik [1, 2]. We extend previous investigations by 
focussing on global aspects of this problem, and use a compactified evolution scheme which includes null infinity on our 
numerical grid. The motivation is twofold: First, we want to simplify the discussion and improve the understanding 
of local versus global issues in critical collapse. In particular, we try to explain questions like: What is the role 
of asymptotic flatness for critical collapse (e.g. the critical solution, the "Choptuon" is self-similar, and thus not 
asymptotically fiat)? How would hypothetical detectors of radiation observe the dynamics close to criticality? How 
can we understand the way null infinity approximates observers at large distances in this simple but nontrivial setup? 
The second motivation is to test numerical algorithms which are based on compactification methods in a situation 
that is very demanding on accuracy. We will argue that at least in the model considered here, global methods do not 
cause a significant penalty in accuracy, but simplify the interpretation of certain results. 

Critical phenomena in gravitational collapse have been originally discovered in the pioneering numerical investiga- 
tions of Choptuik [1, 2]. He studied a massless scalar field coupled to gravity with sophisticated numerical techniques 
that allowed him to analyze the transition in the space of initial data between dispersion to infinity and the formation 
of a black hole. It turned out, that black holes of arbitrarily small mass can be created, and that the critical solutions 
approach a discretely self-similar solution, called the "Choptuon" . Both the Choptuon and the scaling law for the 
black hole mass are universal for arbitrary families of initial data. Hamade and Stewart [3] have found numerical 
evidence, that the critical solution contains a naked singularity which can be seen at future null infinity. Similar 
critical solutions - exhibiting (continuous or discrete) self-similarity - have also been found for several other types of 
matter fields, and have been constructed directly in several cases [4, 5, 6, 7, 8]. The problem has also been studied 
extensively from an analytic point of view by Christodoulou [9, 10, 11, 12], in particular he could show that the space 
of regular initial data that lead to naked singularities has measure zero [13]. 

In the current work, we refer to critical collapse phenomena as "critical collapse at the threshold of apparent 
horizon formation" to avoid possible misunderstandings, since critical collapse is essentially a quasilocal phenomenon 
and the standard definition of black holes is based on global concepts (see textbooks like e.g. [14]). Also, this term 
emphasizes the relation of these phenomena to other areas in nonlinear partial differential equations (PDEs), where 
related phenomena occur, but the concept of black holes is absent. 
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Critical behavior of the kind originally found by Choptuik is usually referred to as type II, because of its formal 
correspondence with type II phase transitions of statistical physics. A different type of critical solutions at the thresh- 
old of black hole formation, corresponding to type I phase transitions, is provided by unstable static configurations - 
like those found by Bartnik and McKinnon [15]. 

Linear perturbation calculations of such critical solutions revealed exactly one unstable mode, which confirmed 
their interpretation as intermediate attractors in the language of dynamical systems. Critical phenomena in general 
relativity are reviewed in [16, 17], including discussions in terms of phase transitions and renormalization group 
techniques familiar from statistical physics. 

A massless scalar field in spherical symmetry exhibits type II critical collapse (there arc no regular stationary or 
time-periodic solutions). Type II critical solutions have been found to exhibit continuous or discrete self-similarity in 
the past lightcone of the singularity. In our case, the critical solution is known to be discretely self-similar (DSS), and 
has been constructed directly as an eigenvalue problem [4]. 

A spacetime is said to be DSS [18] if it admits a discrete diffeomorphism $a which leaves the metric invariant up 
to a constant scale factor: 



where A is a dimensionless real constant and n G N. 

We choose scalar field critical collapse in spherical symmetry for several reasons: the model is very well studied and 
we can compare with a large amount of previous numerical and analytical results. Furthermore, the model is also very 
demanding: The value of the echoing period in the DSS critical solution is A ~ 3.44, which is quite larger compared 
to many other models. Note that larger values of A make it more difficult to resolve a large number of echos. 

Our numerical method is based on a characteristic initial value problem, i.e. we foliate spacetime by null cones. 
This allows for a very efficient evolution system and simplifies the study of the causal structure of the solutions. In 
spherical symmetry, caustics are restricted to the center of symmetry, so wc do not have to deal with the dynamical 
appearance of caustics, which causes potential problems for characteristic initial value problems in higher dimensions. 
Our numerical approach mixes techniques from previous work of Garfinkle [19] and the Pittsburgh group [20, 21], in 
particular we follow Garfinkle in moving along ingoing null geodesies to utilize gravitational focusing for increasing 
resolution in the region of large curvature. Furthermore, compactification methods are well studied and relatively 
straightforward to implement in characteristic codes. 

An important aspect of our compactified characteristic evolution scheme is that at late times our null slices asymp- 
totically approach the event horizon, see Fig. 1. Essentially this is because our coordinates can not penetrate a 
dynamical horizon [22, 23, 24, 25] (they become singular at a marginally trapped surface, e.g. at an apparent hori- 
zon), which is spacelike if any matter or radiation falls through it and null otherwise [24, 26]. Note that the dynamical 
horizon is contained inside of the event horizon, and the outermost dynamical horizon approaches (or coincides with) 
the event horizon at late times, assuming cosmic censorship holds. This fact makes our approach in some sense 
complementary to previous critical collapse studies, which were not adapted to the asymptotic regime. In this paper 
we focus on those aspects of critical collapse which are associated with global structure, and in particular the phe- 
nomenology seen by asymptotic observers. We will only make brief remarks about other well-studied aspects such as 
mass scaling and universality. 

The paper is organized as follows: in Sec. II we give a brief introduction to scalar field critical collapse in spherical 
symmetry, and discuss our geometric setup, which is based on Bondi-type coordinates and a compactification scheme 
which introduces the Misner-Sharp mass function as an independent evolution variable, which renders our evolution 
system regular at null infinity. Our numerical algorithms are presented in Sec. III. In Sec. IV we give a detailed 
discussion of our numerical results, which include the study of power law tails, and compare the radiation signal at 
null infinity with a heuristic estimate based on self-similar scaling and quasi-normal mode frequencies. Our results 
and conclusions are summarized in Sec. V. 



We assume spherical symmetry with a regular center. We introduce a Bondi coordinate system {u, r, 9, (f} on 
spacetime based upon outgoing null hypersurfaces u = constant, with the line element 



($l)"5 = e2"^, 



(1.3) 



II. THE CONTINUUM PROBLEM 



A. Geometric Setup 



\ r 



du -h 2dr 




(2.1) 
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FIG. 1: A Penrose diagram of a typical collapse spacetime. Shown is our numerical null grid which extends to future null 
infinity J^"*". The grid consists of the null slices u = const and ingoing radial null geodesies v — const. Evolution slows down 
in the vicinity of the future event horizon . We also indicate lines (a) r = const < 2Af/, (b) r = const — 2Mf and (c) 
r = const > 2Mf, where Mf is the final black hole mass. 



and assume that spacetime admits a regular center r = 0. Smothness of the metric requires the metric functions 
/3 and V/r to be smooth functions on the spacetime manifold. Note that r is not a smooth function in this sense. 
Consequently, at fixed retarded time uq, (3 and V behave as 



/3(uo,r)-0(r2), 
F(iio,r) =r + 0(r3), 



(2.2) 



where the gauge has been fixed such that the family of outgoing null cones emanating from the center is parametrized 
by the proper time u at the center. 

Following earlier work [20, 27, 28], we write the curved space wave equation (1.2), □g</) = 0, in terms of a rescaled 
field -0 = as 



□ 



-2/3 



r / ,r r 



= 0, 



(2.3) 



where is the 2 dimensional wave operator in the (u, r) submanifold. This ansatz factors out the known falloff of 
at large distances. The evolution system is completed by two hypersurface equations for the metric functions (which 
follow from the (r, r) and (m, r) components of Einstein's equations, respectively) 



(2.4) 



/3,r = 27rr (0_r) , 
Vr = e^^. 

In spherical symmetry, there exists a well defined notion of quasilocal energy, the Misner-Sharp mass- function [29]: 



m(u,r) = - 



1 e-^^ 

r 



(2.5) 



Note that m/r is a smooth function. The Misner-Sharp mass measures the energy content of a sphere fo radius r and 
reduces the Arnowitt-Deser-Misner (ADM) and Bondi masses in the appropriate limits. 
In adapted coordinates 



= -In- 



{u* - u)C{t) C{t)u* ' 



(2.6) 
(2.7) 
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where u* is a real number which denotes the accumulation time of DSS and ({t + = C{t)j we have 

r{T + nA,z) = r{T,z), (2.8) 
where /* denotes (f) or the metric functions (3, y , which are defined in section II A. In addition also satisfies 

($1/2)"'/'= (-1)" -A, (2.9) 
so that fields even in (f), such as (3,^, m, are periodic in r with period A/2. 

B. Compactification 

Compactification is rather simple for characteristic codes, and has been used extensively in the characteristic 

approach to numerical relativity [20, 27, 28]. Before introducing compactification wc want to say a few words about 
asymptotic series expansions. Assuming initial data that are smooth at J^"*", one can expand the massless scalar field 
in powers of 1/r near 0/+ 

</,(u,r) = ^ + ^+0(r-3). (2.10) 
r 

The coefficient Cnp of the 1 /r^- term in the expansion is a Newman-Penrose constant [30] of the scalar field. Inserting 
the expansion (2.10) into the hypersurface equations (2.4) yields 

p{u, r) = H{u) - + 0{r-% (2.11) 

and 

V{u, r) = e2^(") - 2M{u) + + 0{r-^), (2.12) 

where integration constants H{u) and M{u) have been introduced. H{u) indicates redshift since Bondi time Ub is 
related to proper time at the center via the relation 

^=e2«("). (2.13) 

M{u) = lim,.^cc to(m, ^)|M=consi IS thc Bondi mass which is in general not conserved. The Bondi mass-loss equation 
(derived from the uu component of Einstein's equations) states 

^ = -AnN{u^f, (2.14) 

where the news-function is defined as 

fir 

N{u^) = (2.15) 
duB 

In order to write down a compactified evolution system, we introduce a compactified radial coordinate 

x:=^, (2.16) 

so that points at .^"^ arc automatically included in the grid at 3; = 1. Our aim is to regularize our equations at null 
infinity, the coordinate singularity at the regular center can be dealt with in a straightforward manner, see [31]. 

A naive approach of rewriting the hypersurface equations in terms of the a;-coordinate leads to a singular equation 
for the quantity V 

P,^ = 27ra;(l - x){cf>,^f, 

e2/3 (2.17) 



V 



(1-X)2 



FIG. 2: A representative null-paralellogram in the numerical grid made up by two u = const surfaces and two ingoing null 

geodesies v = const. 



To obtain a regular system of evolution equations, we eliminate V by the Misner-Sharp mass-function m{u,r). The 
set of hypersurface equations then becomes: 



dm „ n 

ax 



2(1 -x) 
1 m 



(2.18) 



= 27rx(l - x){^,^f. 



Note that these equations arc completely regular. 

We choose our gridpoints to freely fall along ingoing radial null geodesies x{u) which fulfill 



dx 
du 



^-{l-xfe^^[l-2m^-^^. (2.19) 



Note that the term = m/r does not cause problems because of the smoothness of the metric at the regular 

center (2.2). In section IIIB we will argue that this choice is crucial to resolve DSS phenomena. 
The vanishing of the outgoing null expansion 

e+ = -6-2/5 (2.20) 

r 

on some 2-sphere r = constant means that this sphere is marginally outer trapped. Since this requires /3 to diverge, 
we cannot penetrate apparent horizons in Bondi coordinates. 

III. NUMERICAL ALGORITHM 
A. Evolution 

We choose to evolve the scalar field (p by the diamond sc;lienie outward marching algorithm [21] which can be 
obtained by integrating (2.3) over a null parallelogram E (see figure 2): 

jnni^ = jjudr{Y.) t (3.1) 

Since any 2 dimensional metric is conformally flat and the conformal weights of D/j and the surface area element 
d^x\/—h cancel, the integral /j, □/jV' is equal to the flat space result —2 dudvtp^uv and we obtain 
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FIG. 3: The convergence of the error diagnostic Euur with increasing grid resolution for two near-critical evolutions. Evolution 
(1) uses 5000 gridpoints and p = p*[5000] + 10"^° and evolution (2) uses 10000 gridpoints and p = p*[10000] + 10"^°. The data 
displayed as dots and circles have been sampled. 



In terms of the ^-coordinate this scheme becomes 

r 1 — X 

ipN = ipw + — ips — / dudxe^^mip — —. (3.3) 

The gauge and regularity conditions for the metric functions and the massless scalar field at the origin of spherical 
symmetry become: 

/3(u,x) = Oix^), 

m{u,x)^0{x^), (3.4) 
ip{u,x) = 0{x). 

Our code is based on the "DICE" (Diamond Integral Characteristic Evolution) code, which has been documented 
in [31] (there particular emphasis is given to detailed convergence tests). 

The code is globally second order accurate. Figure 3 shows a convergence test for near-critical evolutions. We 
want to emphasize that the critical value p* of the initial data parameter depends on the grid resolution. This fact is 
essential when doing convergence tests for near critical evolutions, as has been discussed in our previous paper [31]. 

Timesteps are chosen adaptively via a condition 

{V/r)Au < CAr, (3.5) 

where C is a constant, on the order of unity. This restriction on Au is most severe at For various analysis 

purposes, gridfunctions can be evaluated at a: = const locations by determining the nearest gridpoint to the specified 
X — const value and using cubic spline interpolation in a 5-gridpoint interval. 

To monitor the accuracy of the code during runs we use components of the Einstein equations which are auto- 
matically satisfied if the evolution equations hold, such as the following linear combination of the (u, u) and (u, r) 
components 

-Euur = [Guu - 87rT„„) - r^{V/r) (g„, - SirTur) ■ (3.6) 
This can be rewritten in the following form 

Euur = 2e^'^m + Sir 

where / = Since this expression is a linear combination of tensor components, we use a suitably normalized 

quantity E„„r — j^^"'' — , where is the sum of the running maxima (maxi<j fi over a, u = cons<-slice) of the 



2m 

X 



(1-X)> V,x(l-X) 



(3.7) 
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FIG. 4: This figure shows the focussing of ingoing null-geodesics by gravity in the late stages of a slightly supercritical evolution. 
The discretely self-similar dynamics causes the density of the geodesies to increase in a periodic manner. 



absolute values of the individual terms of Euur- We must have Euur 'C 1 for our finite diffcrcnec solution to be a 
good approximation to the continuum solution. Provided that there are no unexpected cancellations between the two 
terms in equation (3.6), Euur will be a reliable measure for the accuracy of our code. A convergence test for this 
quantity is shown in figure 3. 



B. Mesh Reflnment 



Hamade and Stewart [3] have implemented full Berger-Oliger mesh refinement in double null coordinates (without 

compactification) to achieve sufficient resolution to study critical collapse. Garfinke [19] has shown, that this is not 
really necessary - here we follow his approach to increase resolution: Most importantly, we choose our gridpoints 
to follow ingoing radial nuUgeodesics. This leads to a rapid loss of gridpoints in the early phase of collapse, but to 
an accumulation of gridpoints in the region of strong curvature for the late stages of critical collapse (see figure 4). 
Furthermore, when half of the gridpoints have reached the origin we refine the grid and thus obtain a very simple but 
effective form of mesh rc^fiiicniicmt which is a crucial ingredient in the; calculation of critical collapse spacctimcs. In 
previous work [31] we have also tuned our outermost gridpoint to be located just outside of the self-similar horizon 
(SSH). Here we choose to go out all the way to null infinity. The most effective approach in this situation would be 
to just refine the region inside the SSH when half of the gridpoints in this region have reached the origin. While this 
is straightforward to implement, we found the penalty on the resolution that the original condition (the loss of half 
of all gridpoints) causes to be acceptable for the results presented here. 



IV. RESULTS 



A. Identification of critical behavior 



We consider 1-parameter families of initial data cf) = (j)p (uq , x) , such that for small values of p we have dispersion, 
while for large values of p we have black hole formation. It has been found numerically [1, 2] that, for any initial 
data family considered, the evolution of near-critical data approach a universal DSS critical solution. For the present 
paper we already assume universality and restrict ourselves to Gaussian-like initial data 



(f>{uo,x) =^r(a;)^exp — ( 



,r{x) 



Toy 



(4.1) 



where r{x) = This choice makes it easy to compare compactified to uncompactified evolutions using the same 

initial data. All results presented here use ro = 0.7 and a = 0.3 and a radial resolution of 10000 gridpoints. The 
criticality parameter p is identified with the amplitude A. Fig. 5 shows ]V'(t, a;)| for a near-critical evolution with 
amplitude A = 0.144930343980315. 
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FIG. 5: This figure shows a surface plot of I'lpir, x)\ for a near-critical evolution. When the initial Gaussian reaches the origin, 

it is "instantly" (in retarded time u) radiated to future null infinity .J^'^ (located at a:: = 1) by interfering nonlinoarly with the 
field that has not yet reached the origin. Once the evolution has come close to the critical solution, the matter field ip{u, x) = 4>r 
decays exponentially; further self-similar features are thus not visible in this plot. 

As has become common practice, we find near-critical data through a bisection search in p. This procedure yields in 
particular a numerical approximation to the critical value p = which defines the threshold of singularity formation. 

In the bisection search, a number of criteria are possible to distinguish dispersion from collapse. Their equivalence 
can be checked a posteriori when comparing the final result, i.e. subcritical and supercritical solutions close to 
criticality. A typical criterion is to monitor the ratio 2m/ r, where 2m /r = 1 signifies the presence of an apparent 
horizon. This c;ritcrion has been used successfully also in combination with slic;ing conditions, that do not penetrate 
apparent horizons - as is the case in our approach. Numerically, we have used the condition 2rn/r > 0.995 as the 
threshold for apparent horizon formation and for estimating the black hole mass. Remarkably, in practice it turns 
out, that 2m/ r > 0.6 is a sufficient criterion to mark a scalar field evolution as supercritical and thus is useful to 
speed up bisection searches. For practical and historical reasons, this is the approach we have adopted for our code. 

A number of other options come to mind, in particular in our context of evolving out to null infinity, one could e.g. 
monitor the redshift or Bondi mass. In a dispersion evolution, the redshift will decay to zero, while it will approach 
infinity when a black hole forms. Similarly, the Bondi mass will decay to zero with a characteristic tail behavior, as 
is shown in Fig. 6, in a dispersion evolution, and will asymptote to a (positive) constant when the field does not 
disperse. 

In the course of a near-critical evolution, remnants of the self-similar dynamics which occur locally, inside the SSH, 
are radiated to future null infinty. Remarkably, we observe that the imprints of DSS behavior are still present in 
asymptotic quantities such as the Bondi mass and the news function (see figures 7 and 8). 

In addition, in numerical evolutions of supercritical data, the black hole mass has been found to exhibit a universal 
scaling law (see [5, 32, 33]): 

InrriBH = 7ln(p - p*) + *(ln(p - p*)) + const, (4.2) 

where 7 w 0.373 and the function ^ is periodic with period 5A/7 in ln(p — p*). One needs to ensure to measure 
the mass of the outermost component of the horizon in order to obtain the fine-structure of the scaling law shown in 
figure 9. 

This scaling law has been derived by arguments based on linear perturbation theory around the critical solution 
and dimensional analysis however without using a precise definition of the black hole mass [5, 32]. A typical approach 
in numerical simulations based on coordinates that do not penetrate apparent horizons, as far as we are aware, is to 
follow a peak in 2m,/ r until this quantity almost reaches unity, at which point the simiilation is usually slowed down 
by a Courant-Friedrichs-Levy (CFL)-type condition, and then read off the approximate horizon mass at this point. 
We follow this heuristic approach, and are able to reproduce the mass scaling and fine structure. 

At least conceptually, this approach is problematic without quantifying how much mass-energy still remains outside 
of the peak in 2m /r where the horizon-mass is read off. From the perturbation theory argument, one expects scaling 
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ncarcritical 
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FIG. 6: We compare the decay of the Bondi mass in supercritical, near-critical and subcritical evolutions. In the subcritical 
case, the Bondi mass is found to decay for late times with a power-law exponent of approximately —5. 
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FIG. 7: This figure plots the Bondi mass ms against both ub and the adapted time tb for a barely supercritical evolution 
with final black hole mass M/ « 5 X 10"^ The Bondi mass rriB and the mass at the past SSH, mssn, are found to decrease 
exponentially in tb (with an overlayed TB-periodic oscillation with period A/2), once the evolution has sufficiently approached 
the critical solution near the center of spherical symmetry. We also show tu-bxt, the energy present outside of the SSH. 

behavior only for quantities within the self-similar region of spacetime, which does not necessarily extend beyond the 
SSH. We will return to this point in the discussion section V. 

Note that for a near-critical solution, the approximate value of the accumulation time naturally defines an approx- 
imate location (i.e. advanced time) of the SSH. 

For observers at null infinity, a natural time coordinate is Bondi time Ub as defined in eq. (2.13) - Bondi time can 
be identified with the proper time of timelike observers at large distances - see the discussion in Sec. V. In order to 
analyze critical phenomena, we define a time coordinate which is suitably adapted to self-similar critical collapse and 
set 



In- 



(4.3) 



The adapted time coordinate Tb can be defined for spacetimes which are close to the critical solution inside of the 
SSH, so that the value of the accumulation time u* can be determined by a fit to periodicity in Tb- We have used a 
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FIG. 8: We show the news function N(u), first as a function of the natural time coordinate ub of an asymptotic observer and 
also as a function of a suitably adapted time tb = — In "^^ "° where N{tb) is periodic with period A ~ 3.44 after the spacetime 
has come close to the critical solution. Even if the constant is not known, it can be determined by a fit to periodicity in tb. 
Thus, it is possible to observe DSS at J^'^ and to extract the critical exponent A. 




FIG. 9: The fine-structure in rriBH after subtracting a linear fit. The measured period 4.6 is close to the value predicted by 
perturbation theory |A/7 « 4.61. 



fit to periodicity of the news function - N{tb) is periodic with period A ~ 3.44 after the spacetime has come close to 
the critical solution (see figure 8). 

Note that Tb is only an approximate adapted coordinate since it depends on the relation between Bondi and central 
time (2.13). In order to gain some insight into the behavior of Tb(t) consider the simpler case of a continuously 
self-similar (CSS) collapse. We assume that (3 changes only little outside the past SSH, i.e. H{u) w (3ssh{u)- Then, 
since /3ssh(u) = constant in the CSS regime, we find by integrating (2.13) that 

Ub ~ Cu, (4.4) 

where C = e^'^ssn ^j^^^ have chosen the initial condition Ub(0) = 0. Furthermore, it follows from the definition of 
the adapted time coordinates, equations (2.6) and (4.3), that 

Tb « T, (4.5) 

in the CSS regime. Numerically, the deviations from this relation for the DSS case turn out to be quite small. 
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FIG. 10: A conformal diagram of a critical collapse spacetime. In the backwards light cone of the accumulation point the 
dynamics are close to the DSS critical solution. The lightlike boundary is also called the past self-similarity horizon (SSH). 
Depending on whether the initial data are sub- or supercritical, the spacetime will for late times be Minkowski or Schwarzschild. 



DSS behavior in the Bondi Mass and the News Function 



The following argument suggests that the news function is approximately periodic in Tb, as shown in figure 8. 
Assume effective DSS data for the scalar field and the metric functions on the past self-similarity horizon (SSH) (see 
figure 10). Moreover, we assume that the contribution of the right hand side integral in the wave equation for the 
scalar field (3.1) can be neglected outside of the SSH, such that the DSS data on the SSH are linearly propagated to 
without backscattering, i.e. 

lim V(w,r) w Vssh(m)- (4.6) 

r — ^oo 

Furthermore, assume that changes in /3 outside of the SSH are small, so that /3ssh(w) « H{u). It then follows that 

dc{uji) di>ssn{u) 



du 

d 

d^ 



du-B 

(j)sSH{T)C{T)e~'^U* 



-2H 



(4.7) 



_ g-2/3sSH 



— (</.ssH(r)C(T))-0ssH(r)C(T) 



which shows that the news function N{ub) is approximately periodic in r (and in Tb if equation (4.5) holds) with 
period A and satisfies N{t + nA/2) = {-!)"■ N{t). 

In order to determine the behavior of the Bondi mass, we can then rewrite the Bondi mass-loss equation (2.14) 



druB 

dTn 



duB 

dTn 



-47r<e-^^iV^(rB). 



Since N'^{tb) is A/2-periodic in Tq, itIb then takes the following form 

TOB(rB) ~ e"'^='/sSH(TB), 



(4.6 



(4.9) 



where /ssh(tb) is A/2-periodic in Tb- 

This behavior mimics the behavior of the mass-function: We rewrite the mass-function in adapted coordinates 
{t,z), using (2.7), 



V 



{z,r)e 



-2/3(2,r) 



(4.10) 
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and evaluate it at the SSH, which, by a judicious choice of the periodic function C(t) can be chosen to be at = 1 
(since the past SSH is a null surface, one needs to ensure that becomes null at z = 1). We obtain 

mssH(r)=e-V;sH(r), (4.11) 

where /s*sh(''') is periodic with period A/2. 

The periodicity and exponential decay of TOb and rrissn are confirmed by our numerical calculations (see fig. 7) . Note 
that the Bondi mass levels off at roughly 10~^ of the mass in the initial configuration, whereas the mass contained 
within the backwards light cone TOssh continues to scale according to the prediction of critical collapse evolution, 
equation (4.9). Since the difference, i.e. ttjext = — mssa, is almost zero at the initial slice, we conjecture that 
rriEXT is the energy due to backscattering in a critical evolution. 

Note that thext oscillates with r, as can be seen by integrating 

dm „ 2m^ / /. 

- = 2..^(l--)(^j , (4.12) 

from the SSH to J^+. The biggest contribution to ttt-ext will come from the vicinity of the past SSH. 

We observe that TOext contained in a slice close to horizon formation is almost constant for different near critical 
evolutions (down to the numerical limit of fine tuning). Therefore, if this TObxt eventually falls through the horizon, 
then the resulting black hole will have a tiny but finite Bondi mass, no matter how fine tuned the data arc . In the 
very late stages of the evolution, the growing redshift, /3 ^ 00 effectively halts the numerical evolution, while the 
error norm max|Euur(TB > 10)| approaches 10~^. 



C. Quasinormal Modes 

In gravity, quasinormal modes (see [34] for a review) are excitations of a black hole (or a star) satisfying radiation 
boundary conditions. These excitations are, in general, obtained from linear perturbations off a fixed backgroimd, 
together with their associated (complex) eigenvalues. Thus, in a highly dynamical setting, such as in critical collapse 
evolutions, one would not expect to see (identify) quasinormal modes. However, it turns out that for our setting, the 
least damped spherically symmetric mode for scalar perturbations of a Schwarzschild black hole plays a relevant role. 

Perturbation theory [35] gives the following value for the half-period 

^ « 28 A3Mbg. (4.13) 

This mode has previously been detected in supercritical evolutions (far away from criticality) for a self-gravitating 
massless scalar field by Gundlach et al. [36]. 

In the following, we analyze radiation signals for near-critical evolutions, where the notion of a fixed background 
mass does no longer apply. We find that the monopole moment of the scalar field c(wb) shows a damped oscillation 
with exponentially decreasing frequency (sec figure 11). Moreover, the sizes of the half-periods measured from one 
extremum to the next in c(mb) roughly agree with the half-periods obtained from the least damped quasi- normal 
mode (QNM) of a Schwarzschild black hole with a strongly changing "background" mass Mf)g{u^) as shown in figure 
12. Mbgiu^) is obtained by evaluating mji{u^) at the mean value between the extrema of c(ub) (which are inflection 
points of TOb(mb)). 



D. Power-law Tails 



As has been established by Price [37], perturbation fields outside of Schwarzschild black holes die off with an 
inverse power-law tail at late times. In contrast to quasinormal modes this behavior does not depend on the details 
of the collapse process, but only on the asymptotic falloff of the effective potential; (i.e. in a curved spacetime, 
wave propagation is not confined to the light cones, rather waves spread inside the light cones, due to scattering off 
spacetime curvature.) Therefore, tail phenomena can be observed independently of the endstate of the evolution. 
For Gaussian initial data the Newman-Penrosc constant [38] vanishes and for late times perturbation theory [36, 39] 
predicts that the field falls off as oc near J'^ and (f) oc near timelike infinity i"*". 

To be more precise [40, 41], consider a distant static observer at r = const and Bondi time Amb elapsed since 
the "main pulse" of radiation has reached the observer (the duration of the main pulse has been assumed negligible 
in [40, 41]). Null infinity is then found to be approximated by the region where Amb <C r within the context of a 
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FIG. 11: This figure shows the scalar field monopole moment c{ub) for a near-critical (but barely supercritical) evolution. The 
half-periods measured from one extremum to the next roughly agree with the prediction of perturbation theory shown in figure 
12. 




FIG. 12: The exponential decay of the QNM half-periods predicted by perturbation theory is shown with annotated values at 
the midpoints between the points of infiection for the same near-critical evolution as in figure 11. 



perturbative analysis of tail behavior [40, 42, 43]. This regime has been termed the "astrophysical zone" by Leaver 
[40, 42, 43]. 

In addition, the convergence of the perturbation expansion in [40, 41] requires Aug ^ M, where M is the mass 
of the background. In our case, the mass which gives rise to the effective potential is bounded from above by the 
mass of the initial (ingoing) pulse. Mi. Therefore, we demand that Aug ^ Mi. Numerically, we observe tails only for 
Aub > lO^Mj. 

Closeness to timelike infinity, on the other hand, demands Aub S> |r*|, where r* is the usual "tortoise" coordinate 
n = r + 2Af ln(^^^). In figure 13 we show power-law exponents determined by fits of ijj at different r = const 
curves over a series of time intervals for a subcritical evolution. As described in section III we use spline interpolation 
to calculate ip{x — const). The exponents have been determined by fitting the field tp at x = constant against a power 
of Mb in 5 distinct time intervals of the evolution. The domains of validity of the exponents predicted by perturbation 
theory, —2 near and —3 near can be observed here. It is clear that the outermost gridpoints in this evolution 
(using 10000 gridpoints) are indeed located in the "astrophysical zone" since r{x = 0.9995) ~ 2000 3> Aub7 where 
Aub is the Bondi time elapsed since the main pulse of radiation has reached the observer at about Ub ~ 3 (see figure 
6). On the other hand, closeness to timelike infinity i"*" demands that Aug ^ r*. Note that « r for r ^ 2Mi, 
where Mi ^ 0.06 is the initial Bondi mass. It is also apparent that for Aub « r the observers are in between the two 
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FIG. 13: This figure shows power-law exponents for a subcritical evolution and illustrates the domains of validity of the 
predictions of perturbation theory for the two zones: —2 near J^"*" and —3 near i"*". 

zones and the power-law exponents seem to ehangc! smoothly. 

Figure 6 displays the power-law decay, with exponent —5, of the Bondi mass for the same subcritical evolution. 
This behavior can be explained by integrating the Bondi mass-loss equation (2.14) with c oc u^^ in the regime of 
power-law tails. 

Note that in situations involving realistic sources and realistic detectors, power-law tails only play a minor (if any) 
role as an actual signal, but the power-law tail results still can provide hints on how to resolve the very important 
question of whether null infinity is a useful idealisation for gravitational wave detectors. Accordingly, we suggest to 
generalize the term "astrophysical zone" to more general, non-perturbative situations, by reinterpreting - in a very 
loose sense - Aub as a suitably chosen large time scale characteristic of the source. The physical idea is that the 
distance from observers of astrophysical phenomena, e.g. gravitational wave detectors, to the radiation sources is 
very large compared to the time during which substantial radiation from the source can be observed. This is at least 
expected for sources where general relativity is important, as opposed to problems for which a (Post-) Newtonian 
approach and the quadrupole formula are sufficient. An example would be a binary black hole merger in another 
galaxy, which might have a characteristic dynamical time scale of a fraction of a second, and which might be observed 
during several thousand cycles, including a portion of the QNM ringdown. 

V. DISCUSSION 

In this work we have presented numerical constructions of portions of near-critical spherically symmetric spacetimes 

that extend up to future null infinity and asymptote to the event horizon. The simulations are based on a compactified 
approach, where the equations have been regularized in a neighborhood of null infinity by introducing the mass function 
as an additional independent variable. The resolution necessary to resolve critical phenomena is gained through letting 
our gridpoints fall along ingoing null geodesies, following Garfinkle [19]. The grid is repopulated with grid points, 
which are filled with values through interpolations when half of the gridpoints have reached the origin. This is not 
optimal, but ensures sufficient resolution for our current purposes. 

We reproduce standard features of near-critical solutions, such as "echoing" and mass scaling with fine structure. In 
addition we extract radiation at null infinity by computing the news function and find a signal with rapidly increasing 
frequency as measured in Bondi time, which is the natural time coordinate associated with far away observers. In 
order to simplify the analysis of near-critical spacetimes, we have defined an approximate adapted time in terms 
of Bondi time, in analogy to the standard time coordinate which is adapted to the DSS discrete diffeomorphism 
(1.3). This coordinate can be used by asymptotic observers to render the signal from a near-critical collapse (almost) 
periodic. The fact that such a definition actually works out, and makes DSS periodicity manifest in quantities defined 
at null infinity is a non-trivial result, which we explain in Sec. IV A. 

Note also that the amplitude of the news function, shown in Fig. 8 stays fairly constant after the initial transient. 
This feature is clearly universal, as long as the radiation signal is dominated by the DSS collapse, since the system 
can be approximated in the DSS region by a perturbation of the critical solution, and, according to our results, 
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the radiation signal is dominated by the DSS structure of the critical solution. Thus, in section IV B, we neglect 
contributions from scalar field which is far outside the DSS region and does not contribute to the critical collapse 
dynamics. Consequently, in this scenario the essential free parameters determining the radiation signal from an actual 
near critical solution are the number of cycles the solution spends in the neighborhood of the critical solution and 
the length scale at which the solution comes close to the critical solution. The robustness of this scenario, i.e. what 
happens if the initial data are such that there is significant mass outside of the SSH, is outside of the scope of this 
paper and an issue for future research. 

Perhaps the most surprising feature of the radiation signal has emerged from our investigation of QNM's, which 
has been motivated by [36], where the first quasinormal mode is found in collapse evolutions and the question is 
posed as to how QN ringing would change close to criticality. We find that even in very close-to-critical evolutions 
there is a correlation of the radiation signal with the period of the first quasinormal mode, determined from the 
time-dependent value of the Bondi mass, as disussed in Sec. IV C, Figs. 11 and 12. This correlation between 
the radiation of the highly dynamical near-critical solution and the quasinormal mode, that is defined in terms of 
perturbations of a static spacetime, certainly deserves further investigation. This surprising feature might even turn 
out to be a key toward understanding the phenomenon of DSS behaviour in near-critical spacetimes. Our results seem 
to suggest that the effective curvature potential for a DSS self-gravitating field acts as a quasi-stationary background 
for scattering processes which can be approximately described by quasinormal modes of a 1-parameter family of 
Schwarzschild black holes with exponentially decreasing mass. The question of the applicability of QNM-motivated 
estimates is quite relevant for numerical relativity, e.g. when extracting wave forms from binary mergers. In the very 
different context of quasinormal modes of Schwarzschild- AdS black holes, Horowitz and Hubeny [44] have pointed 
out an agreement of the numerical values of the critical exponent and the imaginary part of quasinormal frequencies, 
and have speculated on a connection between critical collapse and quasinormal modes. While Horowitz and Hubeny 
assume pure coincidence as the most likely explanation, it seems worth to keep in mind in future work concerning 
critical collapse or quasinormal modes. 

For subcritical initial data we can evolve for very long times, and thus are able to observe power-law tail behavior 
as shown in Figs. 13 and 6. Analytical calculations predict different falloff rates for radiation along null infinity and 
along timelike lines, and the natural question arises, which falloff rate would be seen by a hypothetical observer (in a 
realistic case, observation of power-law tails would require an extremely large signal-to-noisc; ratio). Accordingly, our 
results depicted in Fig. 13 show how the rates at finite but large radius correspond to the value for null infinity for 
a while before they approach the expected late-time value for finite radius. The interpretation of this phenomenon 
is suggested by perturbative work of [40, 41], where different tail falloff rates are computed. There, different regions 
of spacetime are identified, where certain approximations hold. Within the perturbative regime, results obtained for 
null infinity are valid for what has been termed the "astrophysical zone" by Leaver [40, 42, 43], and which is defined 
as the region where Aub ^ r. (Note, however, that Mb ^ Af must also be satisfied). The physical idea is that the 
distance from observers of astrophysical phenomena, e.g. gravitational wave detectors, to the radiation sources is very 
large compared to the time during which the source radiates at an observable rate. 

We argue that our results illustrate that the relevant falloff from the point of view of an astrophysical observer 
is the falloff rate at null infinity, in accordance with the prediction from perturbation theory. We believe that this 
is a nice model calculation that exemplifies that null infinity is indeed a useful approximation for observers at large 
distance from the source, in the sense that such observers are located in the "astrophysical zone" . 

Along the same lines, we would like to point out that by taking appropriate limits in a conformally compactified 
manifold, worldlines of increasingly distant geodesic observers converge to null geodesic generators of future null 
infinity and proper time converges to Bondi time [45]. Note also that a naive correspondence between observers at 
large distance and spatial infinity would be problematic, e.g. compactification at spatial infinity leads to "piling up" of 
waves, whereas at null infinity this effect does not appear - waves leave the physical spacetime through the boundary 
at null infinity. 

Under practical circumstances, e.g. computing the signal from a source of radiation which is locted at a cosmological 
distance from the detector, null infinity more realistically corresponds to an observer that is sufficiently far away from 
the source to treat the radiation linearly, but not so far away that cosmological effects have to be taken into account. 
We are however not aware of a discussion where this sloppy picture has been made more precise. 

When looking at critical collapse from a global spacetime perspective, as we have done here, one is confronted with 
some issues concerning the mass-scaling, that we would like to comment on briefly: When asking the question of 
whether infinitesimally small black holes can be formed - the question which triggered the original work on critical 
collapse - it could be phrased in two slightly different ways: (i) can we form black-hole solutions where the final-state 
black hole has arbitrarily small mass, or (ii) can we form arbitrarily small apparent horizons. The question which has 
been aswercd in the affirmative by critical collapse research is the second one. The first one still seems open. Our 
results for near-supercritical evolutions seem to indicate that the final black hole is significantly larger than the mass 
leading to scaling. We conjecture that due to backscattering of the outgoing radiation one cannot form arbitrarily 
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small black holes, no matter how fine tuned the data are. 

Finally, we want to emphasize the obvious fact that the failure of the critical solution to be asymptotically flat is 
perfectly consistent with its role in the dynamics of a localized object emitting radiation to null infinity. First, note 
that for the dynamics of critical collapse, only a small region of spacetime is relevant. The scale- invariance of the 
DSS solution is compatible with asymptotic flatness in the following way: a near critical solution, which may or may 
not be asymptotically flat, comes close to the critical solution at a length scale which depends on the initial data, 
follows the critical solution for a number of cycles, and then collapses or disperses. The (discrete) self-similarity of 
the critical solution is on the one hand responsible for the failure of the critical solution to be asymptotically flat but 
actually makes it possible to keep a free length scale in the problem of massless scalar field collapse, as is expected to 
allow for asymptotically flat near critical solutions. 
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